[The Polymorphism Analysis of HLA Class II Alleles Based on Next-Generation Sequencing and Prevention Strategy for Allele Dropout].

OBJECTIVE: To investigate the accuracy of next-generation sequencing technology (NGS) in detecting the polymorphisms of HLA-DRB1, DQB1, DQA1, DRB3, DRB4, DRB5, DPA1 and DPB1 alleles in randomly-selected unrelated healthy individuals from Shenzhen Han population, investigate the potential reason for HLA-DRB1 allele dropout in routine NGS, and establish an internal quality control system. METHODS: NGS-based HLA class II genotyping was performed on 1 012 samples using the MiSeqDxTM platform. The suspected missed alleles indicated by the quality control software and HLA-DRB1 homozygotes were confirmed by PCR-SSOP or PCR-SBT methods. RESULTS: A total of 139 alleles were detected, including HLA-DRB1(45), DRB3(7), DRB4(5), DRB5(7), DQA1(17), DQB1(21), DPA1(10) and DPB1(27). HLA-DRB1*09:01(17.09%),15:01(10.72%); DRB3*02:02(25.99%),03:01(10.18%); DRB4*01:03(36.46%); DRB5*01:01(15.42%); DQA1*01:02(20.01%),03:02(17.19%); DQB1*03:01(19.47%),03:03(17.98%), 05:02(11.66%), 06:01(10.67%); DPA1*02:02(54.45%), 01:03(31.18%) and DPB1*05:01(39.13%), 02:01(16.90%) alleles were the most common alleles in Shenzhen Han population (frequencies >10%). There was no statistical difference between the gene frequencies of HLA-DRB1 and DQB1 loci in our study. The HLA Common and Well-Documented Alleles in China (CWD2.4) (χ2=12.68, P >0.05). 94 cases of HLA-DRB1 homozygous samples detected by NGS were retested by PCR-SSOP or SBT method, and one case of allele dropout at HLA-DRB1 locus was found. SBT method confirmed that the allele of DRB1*04:03 was missed. The laboratory internal quality control system was established. Two cases of new alleles were detected and named by WHO Nomenclature Committee for Factors of the HLA System. CONCLUSION: The HLA genotyping results based on NGS showed a significantly lower ambiguity rate. The HLA class II alleles exhibit genetic polymorphism in the Han population of unrelated healthy individuals in Shenzhen. The independent method based on NGS in clinical histocompatibility testing has limitations and requires internal quality control strategies to avoid allele-dropout events.

HLA Gene Polymorphisms in Romanian Patients with Chronic Lymphocytic Leukemia.

Materials and Methods: This study included 66 patients with CLL, diagnosed between 2020 and 2022, and 100 healthy controls. HLA class I and class II genes (HLA-A/B/C, HLA-DQA1/DQB1/DPA1/DPB1, and HLA-DRB1/3/4/5) were investigated using next-generation sequencing technology. Results: Several HLA alleles were strongly associated with CLL. The most important finding was that HLA-DRB1∗04:02:01 (p=0.001, OR = 1.05) and HLA-DRB3∗02:01:01 (p=0.009, OR = 1.03) have a predisposing role in CLL development. Moreover, we identified that HLA-A∗24:02:01 0.01 (p=0.01, OR = 0.38), HLA-DQA1∗05:05:01 (p=0.01, OR = 0.56), HLA-DQB1∗03:02:01 (p=0.03, OR = 0.40), and HLA-DRB4∗01:03:01 (p=0.03, OR = 0.54 alleles have protective roles. Correlations between HLA expression and gender showed that women had a higher expression of protective HLA alleles when compared to men. Conclusions: Our data are the first to indicate that in Romanian patients with CLL, the HLA-A∗24:02:01 and HLA-DQA1∗05:05:01 alleles have a protective role against CLL development, whereas HLA-DRB1∗04:02:01 and HLA-DRB3∗02:01:01alleles are positively associated with CLL.

Two novel HLA-DRB3 alleles identified in Brazilian bone marrow donors: HLA-DRB3*01:122 and -DRB3*01:123.

We identified two new HLA-DRB3 alleles in Brazilian individuals using next generation sequencing.

Reduced human leukocyte antigen mismatching is associated with more favourable outcomes after unrelated donor haematopoietic stem cell transplantation.

The patient-donor human leukocyte antigen (HLA) match remains the most important prognostic factor for successful unrelated donor haematopoietic stem cell transplantation (UD-HSCT). This single-centre study comprised 125 adult patients with malignant haematological diseases undergoing their first UD-HSCT. The primary goal of this study was to validate the impact of HLA matching on HSCT outcomes, specifically at the HLA-DPB1 and HLA-DRB3/4/5 loci. A multivariable Cox regression analysis with a backward selection algorithm was employed to assess the associations of selected prognostic factors with outcomes after UD-HSCT. Any HLA locus mismatch was found to be associated with an increased incidence of grade II-IV acute graft versus host disease (aGvHD) at 100 days (p = .031; hazard ratio [HR] 1.935) and 6 months (p = .004; HR 2.284) after HSCT. The results of the following analyses also confirmed the strong impact of HLA-DPB1-only mismatch on the incidence of grade II-IV aGvHD at 100-day (p = .006; HR 2.642) as well as at 6-month (p = .007; HR 2.401) time periods. The HLA-DPB1-only mismatch was also shown to be statistically significantly associated with lower relapse incidence (p = .034; HR 0.333). The impact of the HLA-DRB3/4/5 mismatch on outcomes was inconclusive, though the two and more HLA-DPB1 + DRB3/4/5-only mismatches showed a trend towards worse outcomes than a single mismatch. Based on our findings and those of more comprehensive studies, the extended HLA loci typing of patients and donors is suggested to avoid unexpected HLA mismatches during the UD selection.

Characterization of the novel HLA-DRB3*02:194 allele by sequencing-based typing.

HLA-DRB3*02:194 differs from HLA-DRB3*02:02:01:02 by one nucleotide substitution in codon 78 in exon 2.

Full genomic sequence of the HLA-DRB3*03:46 allele identified by third generation sequencing.

Full-length sequence of HLA-DRB3*03:46 covers the 5'-untranslated region (UTR), all introns and exons and the 3' UTR.

Characterization of the novel HLA-DRB3*02:193 allele by sequencing-based typing.

HLA-DRB3*02:193 differs from DRB3*02:02:01:11 by one nucleotide substitution in exon 1, and intronic changes.

Are the integrin binding motifs within SARS CoV-2 spike protein and MHC class II alleles playing the key role in COVID-19?

The previous studies on the RGD motif (aa403-405) within the SARS CoV-2 spike (S) protein receptor binding domain (RBD) suggest that the RGD motif binding integrin(s) may play an important role in infection of the host cells. We also discussed the possible role of two other integrin binding motifs that are present in S protein: LDI (aa585-587) and ECD (661-663), the motifs used by some other viruses in the course of infection. The MultiFOLD models for protein structure analysis have shown that the ECD motif is clearly accessible in the S protein, whereas the RGD and LDI motifs are partially accessible. Furthermore, the amino acids that are present in Epstein-Barr virus protein (EBV) gp42 playing very important role in binding to the HLA-DRB1 molecule and in the subsequent immune response evasion, are also present in the S protein heptad repeat-2. Our MultiFOLD model analyses have shown that these amino acids are clearly accessible on the surface in each S protein chain as monomers and in the homotrimer complex and bind to HLA-DRB1 ß chain. Therefore, they may have the identical role in SARS CoV-2 immune evasion as in EBV infection. The prediction analyses of the MHC class II binding peptides within the S protein have shown that the RGD motif is present in the core 9-mer peptide IRGDEVRQI within the two HLA-DRB1*03:01 and HLA-DRB3*01.01 strong binding 15-mer peptides suggesting that RGD motif may be the potential immune epitope. Accordingly, infected HLA-DRB1*03:01 or HLA-DRB3*01.01 positive individuals may develop high affinity anti-RGD motif antibodies that react with the RGD motif in the host proteins, like fibrinogen, thrombin or von Willebrand factor, affecting haemostasis or participating in autoimmune disorders.

Full genomic sequence of the HLA-DRB3*02:32 allele by Single Molecule Real-time Sequencing Technology.

Full-length sequence covers the 5'-untranslated region (UTR), all introns and exons and the 3' UTR.

Targeted capture enrichment and sequencing identifies HLA variants associated with the severity of COVID-19.

BACKGROUND: Coronavirus disease 2019 (COVID-19) is currently a global pandemic. The pathogenesis of severe COVID-19 has been widely investigated, but it is still unclear. Human leukocyte antigen (HLA) plays a central role in immune response, and its variants might be related to COVID-19 progression and severity. OBJECTIVE: To investigate the hypothesis that individual HLA variations could alter the course of COVID-19 and might be associated with the severity of COVID-19. METHODS: In this study, we conducted an HLA targeted capture enrichment and sequencing of severe COVID-19 patients matched to mild cases. A total of 16 COVID-19 patients, confirmed by SARS-CoV-2 viral RNA polymerase-chain-reaction (PCR) test and chest computed tomography (CT) scan, were enrolled in this study. The HLA targeted capture enrichment and sequencing were conducted. HLA typing was performed by comparing contigs with IPD-IMGT/HLA Database. RESULTS: In this study, 139 four-digit resolution HLA alleles were acquired. The results showed that HLA-DRB3*01:01 allele was significantly associated with the severity of COVID-19 (odds ratio [OR] = 27.64, 95% confidence interval [CI] = 1.35-560.50, P = 0.0064). And HLA-K*01:01 might be a potential risk factor for COVID-19 severity (OR = 0.11, 95% CI = 0.017-0.66, P = 0.019), but HLA-K*01:02 might be a protective factor (OR = 7.50, 95% CI = 1.48-37.92, P = 0.019). CONCLUSION: Three non-classical HLA alleles, including HLA-DRB3*01:01, HLA-K*01:01, HLA-K*01:02 were identified to be associated with the severity of COVID-19 by comparing mild and severe patients. The current findings would be helpful for exploring the influence of HLA gene polymorphisms on the development and severity of COVID-19.

Co-expression network analysis identifies potential candidate hub genes in severe influenza patients needing invasive mechanical ventilation.

BACKGROUND: Influenza is a contagious disease that affects people of all ages and is linked to considerable mortality during epidemics and occasional outbreaks. Moreover, effective immunological biomarkers are needed for elucidating aetiology and preventing and treating severe influenza. Herein, we aimed to evaluate the key genes linked with the disease severity in influenza patients needing invasive mechanical ventilation (IMV). Three gene microarray data sets (GSE101702, GSE21802, and GSE111368) from blood samples of influenza patients were made available by the Gene Expression Omnibus (GEO) database. The GSE101702 and GSE21802 data sets were combined to create the training set. Hub indicators for IMV patients with severe influenza were determined using differential expression analysis and Weighted correlation network analysis (WGCNA) from the training set. The receiver operating characteristic curve (ROC) was also used to evaluate the hub genes from the test set's diagnostic accuracy. Different immune cells' infiltration levels in the expression profile and their correlation with hub gene markers were examined using single-sample gene set enrichment analysis (ssGSEA). RESULTS: In the present study, we evaluated a total of 447 differential genes. WGCNA identified eight co-expression modules, with the red module having the strongest correlation with IMV patients. Differential genes were combined to obtain 3 hub genes (HLA-DPA1, HLA-DRB3, and CECR1). The identified genes were investigated as potential indicators for patients with severe influenza who required IMV using the least absolute shrinkage and selection operator (LASSO) approach. The ROC showed the diagnostic value of the three hub genes in determining the severity of influenza. Using ssGSEA, it has been revealed that the expression of key genes was negatively correlated with neutrophil activation and positively associated with adaptive cellular immune response. CONCLUSION: We evaluated three novel hub genes that could be linked to the immunopathological mechanism of severe influenza patients who require IMV treatment and could be used as potential biomarkers for severe influenza prevention and treatment.

Association between germ-line HLA and immune-related adverse events.

Background: In recent years, significant progress has been made in immune checkpoint inhibitors (ICIs). However, accompanied by remarkable efficacy, a growing number of immune-related adverse events (irAEs) also arose. The mechanism of irAEs remains unclear. Previous studies indicated a positive association between specific human leukocyte antigen (HLA) variants and irAEs. Therefore, we planned and initiated a large cohort study aiming to uncover the relationship between irAEs and divergent HLA types. Methods: We screened all patients who have been treated in the clinical research ward, Cancer Hospital of the Chinese Academy of Medical Sciences. All participants were diagnosed with malignant tumors with complete AE follow-up data in the original electronic medical records. Sequencing libraries were generated using a customized panel, and four-digit formatted HLA alleles were extracted for further analysis. Association analysis was performed between HLA variants and different irAEs. We introduced two external reference groups and a non-irAE control group within the study cohort to control the type I error. We also explored the relationship between the zygosity of HLA genes, the evolutionary divergence of HLA class I genotype (HED), and irAEs. Results: 530 participants received at least two doses of ICIs. The median follow-up time was 10.3 months. 97% of patients received anti-PD-1/PD-L1 treatment. The occurrence of overall irAEs showed no significant difference between the HLA homozygous group and the HLA heterozygous group. We did not find any significant association between irAEs and HED. We found that some HLA types are associated with irAEs of different organs and detected a significant association between HLA-DRB3*01:01 and thrombocytopenia (OR 3.48 (1.19,9.42), p = 0.011), HLA-DPB1*04:02 and hypokalemia/hyponatremia (OR 3.44 (1.24,9.1), p = 0.009), leukopenia (OR 2.1 (0.92,4.8), p = 0.037), anemia (OR 2.33 (1.0,5.41), p = 0.026), HLA-A*26:01 and bilirubin elevation (OR 2.67 (0.92,8.31), p = 0.037). Conclusions: IrAEs in specific organs and tissues may be associated with certain HLA types, while HLA heterogeneity has no significant influence on the happening of irAEs. More research is needed to explore the role of germline genetic changes in the risk assessment of irAEs.

Characterization of the novel HLA-DRB3*02:179N allele by sequencing-based typing.

HLA-DRB3*02:179N differs from DRB3*02:02:01:02 by one nucleotide substitution in codon 98 in exon 3.

High-resolution HLA class II sequencing of Swedish multiple sclerosis patients.

Multiple sclerosis (MS) is a chronic neurological disease believed to be caused by autoimmune pathogenesis. The aetiology is likely explained by a complex interplay between inherited and environmental factors. Genetic investigations into MS have been conducted for over 50 years, yielding >100 associations to date. Globally, the strongest linkage is with the human leukocyte antigen (HLA) HLA-DRB5*01:01:01-DRB1*15:01:01-DQA1*01:02:01-DQB1*06:02:01 haplotype. Here, high-resolution sequencing of HLA was used to determine the alleles of DRB3, DRB4, DRB5, DRB1, DQA1, DQB1, DPA1 and DPB1 as well as their extended haplotypes and genotypes in 100 Swedish MS patients. Results were compared to 636 population controls. The heterogeneity in HLA associations with MS was demonstrated; among 100 patients, 69 extended HLA-DR-DQ genotypes were found. Three extended HLA-DR-DQ genotypes were found to be correlated to MS; HLA-DRB5*01:01:01-DRB1*15:01:01-DQA1*01:02:01-DQB1*06:02:01 haplotype together with (A) HLA-DRB4*01:01:01//DRB4*01:01:01:01-DRB1*07:01:01-DQA1*02:01//02:01:01-DQB1*02:02:01, (B) HLA-DRBX*null-DRB1*08:01:01-DQA1*04:01:01-DQB1*04:02:01, and (C) HLA-DRB3*01:01:02-DRB1*03:01:01-DQA1*05:01:01-DQB1*02:01:01. At the allelic level, HLA-DRB3*01:01:02 was considered protective against MS. However, when combined with HLA-DRB3*01:01:02-DRB1*03:01:01-DQA1*05:01:01-DQB1*02:01:01, this extended haplotype was considered a predisposing risk factor. This highlights the limitations as included with investigations of single alleles relative to those of extended haplotypes/genotypes. In conclusion, with 69 genotypes presented among 100 patients, high-resolution sequencing was conducted to underscore the wide polymorphisms present among MS patients. Additional studies in larger cohorts will be of importance to define MS among the patient group not associated with HLA-DRB5*01:01:01-DRB1*15:01:01-DQA1*01:02:01-DQB1*06:02:01.

Sib-pair subgroup familial type 1 diabetes mellitus in children in the state of Qatar.

BACKGROUND: Type 1 diabetes is the most common type of diabetes mellitus (DM) in children. It can be sporadic in onset or cluster in families, which comprises parent-offspring and sib-pair subgroups. The risk of developing DM in first-degree relatives of affected individuals is 8-15 fold higher. There is limited data about familial DM from the Gulf region. This study aims to describe the clinical, biochemical and genetic characteristics of sib-pair familial type 1 diabetes in Qatar. METHODS: Every child with DM following up at Sidra Medicine was recruited. Data was collected regarding clinical features, family history, type 1 diabetes autoantibodies and whole genome sequencing was performed. Genetic analysis for MODY genes and HLA association analysis was conducted. RESULTS: 44 families with sib-pair familial diabetes were identified. Of these, 2 families had 4 affected siblings and 5 families had 3 affected siblings. The majority are of Qatari ethnicity and the most common autoantibody was GAD65. The most common age of onset in the proband was 5-9 years while it was 10-14 years in subsequent siblings. The occurrence of DKA & HbA1c levels were lower in the second affected sibling. No relevant MODY gene variants were found. HLA analysis found 15 variants in at least 50% of the subjects. Most common were HLA-F*01*01*01G, HLA- DPA1*01*03*01G, HLA- DRB3*02*02*01G, HLA- E*01*01*01G & DRB4*03*01N. CONCLUSIONS: The prevalence of sib-pair diabetes is 3.64%. The second affected siblings were older. MODY is unlikely and Class I and II HLA genes was present in sib-pair diabetes.

Full genomic sequence of the HLA-DRB3*02:22:01 allele by single molecule real-time sequencing technology.

Full-length sequence covers the 5'-untranslated region (UTR), all introns and exons, and the 3' UTR.

Identification of the novel HLA-DRB3*02:174 allele by next-generation sequencing.

HLA-DRB3*02:174 differs from DRB3*02:02:01 by one nucleotide substitution in codon 95 in exon 3.

Next-generation sequencing and genotype association studies reveal the association of HLA-DRB3*02:02 with delayed hypersensitivity to penicillins.

BACKGROUND: Nonimmediate (delayed)-allergic reactions to penicillins are common and some of them can be life-threatening. The genetic factors influencing these reactions are unknown/poorly known/poorly understood. We assessed the genetic predictors of a delayed penicillin allergy that cover the HLA loci. METHODS: Using next-generation sequencing (NGS), we genotyped the MHC region in 24 patients with delayed hypersensitivity compared with 20 patients with documented immediate hypersensitivity to penicillins recruited in Italy. Subsequently, we analyzed in silico Illumina Immunochip genotyping data that covered the HLA loci in 98 Spanish patients with delayed hypersensitivity and 315 with immediate hypersensitivity compared to 1,308 controls. RESULTS: The two alleles DRB3*02:02:01:02 and DRB3*02:02:01:01 were reported in twenty cases with delayed reactions (83%) and ten cases with immediate reactions (50%), but not in the Allele Frequency Net Database. Bearing at least one of the two alleles increased the risk of delayed reactions compared to immediate reactions, with an OR of 8.88 (95% CI, 3.37-23.32; p < .0001). The haplotype (ACAA) from rs9268835, rs6923504, rs6903608, and rs9268838 genetic variants of the HLA-DRB3 genomic region was significantly associated with an increased risk of delayed hypersensitivity to penicillins (OR, 1.7; 95% CI: 1.06-1.92; p = .001), but not immediate hypersensitivity. CONCLUSION: We showed that the HLA-DRB3 locus is strongly associated with an increased risk of delayed penicillin hypersensitivity, at least in Southwestern Europe. The determination of HLA-DRB3*02:02 alleles in the risk management of severe delayed hypersensitivity to penicillins should be evaluated further in larger population samples of different origins.

HLA-DRB3/4/5 Matching Improves Outcome of Unrelated Hematopoietic Stem Cell Transplantation.

The HLA-DRB3/4/5 loci are closely linked to the HLA-DRB1 gene. Mismatches in these loci occur with a frequency of about 8%-12% in otherwise 10/10 HLA-matched transplant pairs. There is preliminary evidence that these disparities may associate with increased acute graft-versus-host disease (GvHD) rates. The aim of this study was to analyze a large cohort of German patients and their donors for HLA-DRB3/4/5 compatibility and to correlate the HLA-DRB3/4/5 matching status with the outcome of unrelated hematopoietic stem cell transplantation (uHSCT). To this end, 3,410 patients and their respective donors were HLA-DRB3/4/5 and HLA-DPB1 typed by amplicon-based next-generation sequencing (NGS). All patients included received their first allogeneic transplant for malignant hematologic diseases between 2000 and 2014. Mismatches in the antigen recognition domain (ARD) of HLA-DRB3/4/5 genes were correlated with clinical outcome. HLA-DRB3/4/5 incompatibility was seen in 12.5% (n = 296) and 17.8% (n = 185) of the 10/10 and 9/10 HLA-matched cases, respectively. HLA-DRB3/4/5 mismatches in the ARD associated with a worse overall survival (OS), as shown in univariate (5-year OS: 46.1% vs. 39.8%, log-rank p = 0.038) and multivariate analyses [hazard ratio (HR) 1.25, 95% CI 1.02-1.54, p = 0.034] in the otherwise 10/10 HLA-matched subgroup. The worse outcome was mainly driven by a significantly higher non-relapse mortality (HR 1.35, 95% CI 1.05-1.73, p = 0.017). In the 9/10 HLA-matched cases, the effect was not statistically significant. Our study results suggest that mismatches within the ARD of HLA-DRB3/4/5 genes significantly impact the outcome of otherwise fully matched uHSCT and support their consideration upon donor selection in the future.

Identification of the novel HLA-DRB3*02:02:19 allele.

The HLA-DRB3*02:02:19 allele differs from DRB3*02:02:01:02 by a single nucleotide change in exon 2.